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ABSTRACT 

In the present work detailed calculations have been carried out on the rotation of polarization vector 
of an electromagnetic wave due to the presence of gravitational field of a rotating body. This has been 
done using the general expression of Maxwell’s equation in the curved space-time. Considering the 
far field approximation (i.e the impact parameter greater than the Schwarzschild radius and rotation 
parameter), the amount of rotation of polarization vector as a function of impact parameter has been 
obtained for a rotating body (considering Kerr geometry). Present work shows that, the rotation of 
polarization vector can not be observed in case of Schwarzschild geometry. This work also calculates 
the effect, considering prograde and retrograde orbit for the light ray. Although the present work 
demonstrates the effect of rotation of polarization vector but it confirms that, there would be no net 
polarization of electromagnetic wave due to the curved space-time geometry in Kerr field. 
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1. INTRODUCTION 


It was way back in 1958, when Balazs (1958) for the first time had shown that the gravitational field of a rotating 
body can actually rotate the direction of polarization vector of an electromagnetic wave passing in its field. The author 
gave an idea about how the gravitational field effects the orientation of the plane of polarization of the electromagnetic 
wave. The basic concept behind the idea of Balazs (1958), was to derive the four Maxwell’s equation ^ 'ckson (1962); 
Born & Wolf (1999); Stcphani (2004); Landau & Lifshitz (1998) and express the electric displacement vector in terms 
of electric and magnetic field. He tried to figure out whether the gravitational field would have any effect on polarization 
or not? It should be noted here, that the work by BalazsBalazs (1958) was carried out before the establishment of 
Kerr metric, so he described the metric as: 
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where L represents the angular momentum of rotating mass and R represents the radius vector of rotating body. He 
investigated the problem further by considering the rotation of a frame instead of a rotating body. 

In the very next year another investigator Plebanski (1960), had put forwarded his investigation on the propagation 
of electromagnetic wave affected by gravitation. He tried to understand the similar physical phenomena of rotation of 
polarization vector due to the gravitational effect in an isolated physical system or macroscopic medium. To resolve the 
physical problem, he took the geometrical optics approximation to estimate the displacement vector and considered 
the case of very far field approximation. In that work from the change of the energy momentum tensor, he proved that 
energy conservation remains intact for the modified electric displacement vector. He described the electric displacement 
vector as: 


-Di — Ei -\- EikEk + SikigkEli 


Where Di represents the displacement vector. He showed that when space becomes flat at very far distance, then 
Di ~ Ei. This work also concentrated on the effect of the factors £ik and gk on the rotation of polarization vector. 
The lowest order approximation for £ik and gk had been shown in terms of matter tensor (as described in Plebanski 
(I960)) and concluded that the effect is very small. Finally, the work gave the total rotation of polarization vector for 
a rotating frame and the author applied the result for some selected cases. This work clearly established the rotation 
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of the polarization vector during the deflection of light by gravitational field and also estimated the amount of rotation 
of polarization vector.In the year 1977 Pineault & Roeder (1977) estimated the rotation of the polarization vector with 
reference to a distance observer and calculated the image distortion due to the effect of the rotation of polarization 
vector. In 1980, Su & Mallett (1980) gave a comprehensive review of the work done till date, to calculate the amount 
of rotation which the polarization vector of a light ray will undergo when passing close to a massive rotating body (viz. 
a black hole). In their review they discussed the cases, where some authors got non-vanishing results and others got 
vanishing results, with the reasons for such differences. 

In a more recent work in last decade, Sereno (2004) investigated the rotation of polarization vector due to gravitational 
field and gave the approximate result for gravitational Faraday rotation. He took the weak field approximation limit 
and analytically presented the result for the gravitational Faraday rotation. The author claimed the model to be 
suitable for estimation of Faraday like rotation of polarization vector up to third order correction to gravitational 
potential. He also discussed briefly the case for a gravitational lens moving with constant velocity. Prior to that, 
a very systematic treatment on the image formation due to a moving gravitational lens was presented by Paczynski 
(1986) 

As discussed above, it has been already well established that the direction of polarization vector gets rotated due 
to the effect of gravity Carroll (2004); Balazs (1958); Vilenkin (1984); Fernandez & Davis (2011). In general, space- 
time curvature created due to the presence of rotating massive body, can cause the rotation of polarization vector of 
electromagnetic wave or light. The Sec.2 of this work, gives the general expression for the rotation of the polarization 
vector for a linearly polarized light ray considering the far field approximation, which means the impact parameter is 
much greater than Schwarzschild radius and rotation parameter. The expression for the electric displacement vector 
in the present work has been written with a knowledge of Kerr metric (Kerr (1963)), which was not possible by Balazs 
(1958) and Plebanski (1960) , as Kerr metric expression was published later. As discussed by Sereno (2004) the 
model has significance in the micro-lensing phenomena, where the massive bodies act as moving micro-lenses in case 
of galactic gravitational lensing of light ray. However, this situation clearly differs from our case of interest, where 
we studied for a single ray the effect on the direction of polarization vector, by a rotating gravitational body. The 
authors Su & Mallett (1980) developed their technique to calculate the rotation angle for two types of observers: in 
the first case the observer is in a global inertial frame (GIF), ie in an asymptotic flat space and in the second case, the 
observer is a non-inertial one, at rest in a locally non-inertial frame(LNRF). The result obtained by them under the 
first ca,se(viz. GIF) can only be compared to the findings from our present study, as we limit our work to an observer 
in an asymptotic flat space only (i.e. at infinite distance from the massive rotating body). In Sec.5, the degree of 
polarization as may be expected due to this phenomena, has been obtained with a null result. To derive the effect of 
gravitational field, here in the present work we derive the relation between the electric displacement vector and the 
electric field. A light source placed at asymptotically flat space emits electromagnetic wave, which travels through 
curved space-time and is received by an observer at another asymptotically flat space. The effect of the gravitational 
field on the electromagnetic wave has been calculated in the present work from the knowledge of the metric tensor 
representing the curved space-time. For simplicity the emitted light from the source has been considered to be linearly 
(or plane) polarized for our convenience. 

2. GENERAL EXPRESSION FOR THE ROTATION OF POLARIZATION VEGTOR OF 
ELEGTROMAGNETIG WAVE IN CURVED SPACE-TIME 

To derive the expression for desired rotation of polarization vector of electromagnetic wave, the expressions for the 
electromagnetic displacement vector have to be calculated in generalized form, from the co-variant form of Maxwell 
equation. There expression can be written down as Landau & Lifshitz (1998); Stcphani (2004); Misner et al. (1973): 

dx^ dx^ dx^ 

(lb) 

After the general expression of the displacement vector and magnetic induction have been established, the general 
expression for rotation of polarization vector in curved space-time for a known gravitational field can be obtained. 

2.1. Expression for the electromagnetic displacement vector in curved space-time 

After writing the Maxwell equations in co-variant form, we may write down the expressions for electric field and 
displacement vector in the corresponding three dimensional form as below Landau & Lifshitz (1998); Frolov & Shoom 
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Ea = Fqol Bap = Fap ( 2 ) 

where E is electric component, B is magnetic component, D is the electric displacement vector and H is the magnetic 
induction of the electromagnetic wave. The a and /3 index run from 1 to 3 representing the spatial axis of the 
four dimensional space-time geometry. In this present work, D and E are also represented in vectorial form in three 
dimensional geometry. Now from Eqn.(2) as outlined by Landau & Lifshitz (1998); Frolov & Shoom (2011) we may 
re-write, the following expression in a more detailed manner. 
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With the help of the identities of unit tensor Landau & Lifshitz (1998); Weinberg (1972); Frolov & Shoom (2011) in 
three dimensional coordinate system, the Epn.).-*) can be re-written in vector form as follows: 

D = +Hxg (4) 

So from Eqn.(4), one can find out the components of the displacement vector for a given gravitational field and 
subsequently the rotation of the polarization vector can be found out for a linearly polarized electromagnetic light ray. 

2.2. Rotation of polarization vector of an electromagnetic wave in the gravitational field of a rotating object 

2.2.1. Light travelling along eguatorial plane 

Let us assume that a linearly (or plane) polarized light ray is coming from an asymptotic flat space(represented 
by r=—oo) and passing through the gravitational field of a rotation object represented by the metric g^j,^. After that 
the ray is received by an observer at r=c». The displacement vector defined by Eqn.( ) for the three dimensional 
coordinate system, in the present case has to be derived. For our convenience, we define the geometry of the system 
as follows: The rotating gravitational mass is placed at the origin of three dimensional cartesian co-ordinate system, 
with rotation axis defining the Z axis. The incoming light ray travelling along a plane which is perpendicular to Z-axis. 
The direction of incoming light ray and the origin of co-ordinate system together define the X-Y plane. The light ray 
approaches the origin with impact parameter ‘b’. The original direction of light ray marks the direction parallel to X 
axis (as illustrated by Fig.(l)). The Y-axis is now defined as normal to X and Z axis. 

In this geometry light ray is contained, within the equatorial plane of rotating body. The incoming light ray assumed 
to be plane polarized and it will have its electric vector projected into the directions (Y,Z) with components Ey and 
Ez. Thus the position angle of polarization (xeq) can be expressed by the relation: 

Xeq = arctan(^) 


( 5 ) 
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Now as the light propagates through the gravitational field, the effect of gravitation can be equated with the change 
in the refractive index of the medium Balazs (1958); Plebanski (1960); Sen (2010); Roy & Sen (2014). As a result in 
the new material medium, the strength of electric field vectors will be now estimated by the displacement vector Dy 
and Dz. Thus the position angle of the polarization vector will be now estimated by the formula, 

Xeg = arctan(^) (6) 

So the rotation of the polarization vector at each location of space, would be some over the total change from the 
position angle defined by Eqn.(5). Now Eqn.(6) only gives the position angle of the polarization vector at any point 
under the influence of gravitational field. To find out the total rotation of the polarization vector at any point, we 
have to calculate the rate of change of displacement vector at a particular point. Now one can find the values of the 
components of the displacement vector along the Y and Z axes, with the help of Eqn.( ) and Eqn.(5). We note that 
the electromagnetic wave is propagating parallel to the X axis and there will be no component of displacement vector 
or the magnetic induction vector along the direction of propagation (which is X axis). So one can write: 
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Combining with Eqn.(u): 
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arctan(^)=arctanf^f^^^) 
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where x^q is the position angle of polarization vector at a given point, under the influence of gravitational field. To 
find the change of electric displacement vector, we must know, the components of electric vector are Ey and Ez along 
the Y and Z axes at that at point r = —oo. Now, if we assume that the polarization vector makes an angle Xeq with 
Y axis at the point r = —oo, then the components of electric vector will be related by the equation Ey = ^eqEz (say), 
^eq is given by, 


E 

Xeq = arctan(^) = arctan(^eq) 


( 9 ) 


Now the different components of magnetic field (H) at the starting point of incoming ray can be written in terms of 
the electric field by the following relation Hy = n x Ey and Hz = fi x Ez, (page-112 of Landau & Lifshitz (1998)), 
where n is the normal unit vector along the direction of propagation(in this case parallel to X axis). So taking only 
the magnitude of magnetic field from Eqn.(7a) and Eqn.(7b) and taking note of our earlier assumption that Ey = ^Ez, 
we may write. 


Dy = Ez{^+gx} 
Xffoo 

Dz = Ez{—^ - ^eqgx) 
V5oo 


(10a) 
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As a result we may finally write. 
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Eqn.(ll) gives the position angle of the rotated polarization vector of the electromagnetic wave. Now at initial point 
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(r = —00 in our case) Dy = Ey and = E^. The general equation for Dy, at any point r(say) it will be, 

Dy= ^p^ +E,g^{r) (12a) 

y9oo{r) 

D, =- E,^,yg,ir) (12b) 
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From Eqn.(12a) and Eqn.(12b) it is clear that both Dy and are function of position only. 

2.2.2. Light propagation along rotation axis 

If the light ray propagates along the rotation axis i.e. along the defined Z axis, then the will have its electric vector 
projected into the directions(X,Y) with components E^ and Ey. Thus the position angle of polarization (xax) can be 
expressed as: 

Xax = arctan(^) (13) 

Ey 

Now after the light ray propagates out through the gravitational field, the strength of electric field vectors will be 
estimated by the displacement vectors Balazs (1958); Plebanski (1960); Sen (2010); Roy & Sen (201-1) Dx and Dy. 
Thus the position angle of polarization vector will be estimated by the formula, 

Xax = arctan(^) (14) 

JJy 

As it has been shown in previous Sec. 2.2.1 one can find the values of Dx and Dy,with the assumption that Ex = ^axEy 
at any point r(say), by the following relations 


Dx = + Eyg,,{r) (15a) 

V900 [r) 

Dy = - Ey^ax9z{r) (15b) 
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Here again from Eqn.(15a) and Eqn.(15b), it is clear that both Dx and Dy are function of position (r) only. 

3. METRIC ELEMENT FOR LIGHT RAY TRAVELLING ALONG AND PERPENDICULAR TO EQUATORIAL 

PLANE 

In this section the effect of gravitational field of a rotating object on the total amount of rotation of polarization 
vector will be discussed. In our case the electric displacement vector as in Eqn.(IOa) and Eqn.(lOb) (and Eqn.l5a and 
15b) have been defined in cartesian co-ordinate system, so the situation demands that we must express the metric g 
or the corresponding line element in cartesian form. The authors Su & Mallett (1980) had suitably parameterized the 
Kerr line element, to look for symmetry and to find a solution to the problem. They considered two geometries, the 
light ray passing parallel to the rotation axis and ray passing along the equatorial plane. The Kerr field line element 
is generally expressed in spherical co-ordinates as given bellow Kerr (1963); Visscr (2007), 


ds^ = ~ “ asiT? Odcj)]^ + + a^)d(j) — adt]"^ + (16) 

ZA 2-1 

where A = — rgV + a? and E = cos^ 9. Bellow we shall consider the above two cases of geometry and evaluate 

the metric in cartesian form for the two cases separately. 


3.1. Light ray travelling along equatorial plane 
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In this case we have taken the light ray to be contained in the equatorial plane of the system, so sin0 = 1 and 
COS0 = 0. Taking this into account, we have S = Thus the above Eqn.(16) can be rewritten as: 

To 


ds^ = —(1 — 

+ (- 
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So the metric for the Kerr line element will be, 
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Here from Eqn.( 18) the goo and g^, (or goa;) can be easily identified as: 
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The spatial component 7 a/ 3 ,( Landau & Lifshitz (1998)) could be extracted from the Eqn.(17) and written as: 
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For the geometry as explained in Fig.(l) for impact parameter b, we can write (approximately): 

r2 = *2 (22) 

As the amount of deflection is small, we have used y=b for entire trajectory of light. We also note r ^ b. Applying 

the result for sufficiently large r in Eqn.(12a) and Eqn.(12b), we can say that displacement vector(D) is a function of 
X only. So the displacement vector at point (x+Ax) will be, 



Thus the change of Dy is Dy{x + Ax) — 
change would be given as: 
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Dy{x) for the change in position co-ordinate by Ax. Therefore, the rate of 
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We note here Ey and are independent of x. Now after some steps of simplification, the rate of change of electric 
displacement vector along the Y axis can be calculated as: 


dDy 

dx 


= E, 


C_ X _ ^ 1 

\/ 'Jx^ — Vg 


^/x'^ +b'^ — r„ 


— Srgab 


+ 6^)1 


(25) 


It is easy to observe that the function is an even function so the total change of electric displacement vector from 
initial point (r = —oo) to final point (r = oo) will be, 


_oo = 2 / E, 
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From Eqn.(26), we can get the electric vector at the final point (r = +oo), by adding the total change {ADy\'^^_^) 
to the initial value of electric vector and it is given as: 


Dy{x = +oo) 


Ey + AEy 
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Similarly from Eqn.(23b), we can find the expression for electric vector along Z axis and it can be expressed as: 


D^x 
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It is also clear from the definition of the problem that Dy{x = oo) and Dz{x = oo) are nothing but the components 
of electric field as will be recognized at asymptotically flat space (x = oo). So the total rotation of the polarization 



















































vector for the light ray travelling along the equatorial plane would be(combining Eqn. (9,11,27,28)) 
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(29) 


The Eqn. (29) gives the general expression for rotation of polarization vector of a light ray travelling along the equatorial 
plane of a Kerr mass, where Xeq is the initial value of the position angle of polarization (electric vector) when light 
leaves its source at r= -oo. We can clearly see when we substitute a = 0 (representing Schwarzschild field), the value 
of Axeg becomes zero from Eqn. (29). This confirms that, for Schwarzschild geometry there will be no rotation of 
polarization vector. 


A 



Massive body 

Figure 1. Shows schematic diagram of an electromagnetic wave passing close to a gravitational mass with impact factor ‘b’ 

3.2. Light ray travelling parallel to the rotation axis 

Now, we consider that the light ray propagates parallel to the rotation axis i.e. the Z axis according to our definition. 
In this case we could consider that the light ray is in the Y-Z plane and the choice of the X axis becomes arbitrary. 
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Now looking back to Boyer-Lindquist equation; Eqn.(16) one would have dcj) = 0 and the Eqn.(16) will be reduced to: 

,2 a 


ds^ 


^dr^ 
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A , 2 
-dt 


Sin 
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where S and A are previously defined in Eqn.(16). Following the same steps as shown in Sec. 3.1 one could simplify 
Eqn.(30) and it is given by: 
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Now following the same assumption as in Sec. 3.1, that for the large value of r(^ b, impact parameter) we can write: 

2 2 I 1,2 

r = z + b 

But from Eqn.(31) we see that here gzo = = 0, so from Eqn.(14) with the help of Eqn.(15a) and Eqn.(15b) we can 

write that 

Ex 

Xax = a.rcta.n{-^) = Xax (32) 

Ey 

So there will be no change of position of polarization angle (y), when electromagnetic radiation propagates parallel to 
the rotational axis of a rotating mass. 

In the work authored by Balazs (1958) it had been shown that the rotation of of polarization vector is non-zero when 
light ray travels parallel to the rotation axis in Kerr field. Thus, our result here is definitely contradicting and totally 
opposite to what Balazs (1958)reported in 1958. However Balazs (1958) did not calculate the value of rotation of 
polarization vector when light ray is travelling along the equatorial plane. 

4. SOME SAMPLE CALCULATIONS FOR THE ROTATION OF POLARIZATION VECTOR; WHEN LIGHT 

RAY TRAVELS ALONG EQUATORIAL PLANE 

In order to study the variation of Ayeg with respect to Xeq, we take its first order derivative from Eqn.(29). Thus 

dAxe 


we see Axeq becomes extremum when, = 0. Thus we have 
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From Eqn.(.-..-i) we can notice that secxeq can not be zero and solving the other part of equation to zero we have. 
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From Eqn.(34) we can find out the value of Xeq for which the rotation in position angle will be extremum. Taking a 
second derivative we can show that, this value represented by Eqn.(34), gives the condition for maximum. 

4.0.1. Case I: Ey = 0 i.e. ^ = 0 or Xeq=0 

If we consider Ey = 0 then ^ will be zero. Now from Eqn.(29) the total rotation of polarization vector will be zero. 

4.0.2. Case 11: Ey = E^ i.e. ^ = 1 or Xeq=45° 

Considering the both components are equal the initial position angle of polarization vector will be 45°. In this case 
from Eqn.(29) we could write: 

n rga 

Axeq = arctan ■ 


) Vb 


- 1 


(35) 
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Figure 2. Shows a plot of Axeij against Position angle {Xeq = arctan(^)) where Xeg represents the orientation of initial 
polarization vector at r = —oo. The Axeg is the rotation in position angle. The three curves represent three plots with impact 
parameter equal to 1, 1.5, and 2 times the Physical Radius. 


Sun is rotating with an angular speed 1.664 x 10“"^ ° /s (degrees per second) and it has moment of inertia 5.8 x 10^® 
kg. rn?. Accordingly, we can calculate the rotational parameter ‘a’ for Sun as 1.69749 km. Considering the impact 
parameter ‘b’ of Sun as the physical radius of Sun (0.7x 10®km.), from Eqn.(35)we can estimate the value of Axeq 
for Sun as 4.2873 micro-arcsec. One can say, there will be no observable change for the rotation of polarization 
vector in case of Sun. However, when we apply our Eqn.( 35) to other fast rotating compact objects, we find the 
result is quite fascinating. For these compact rotating objects, the light ray passes closer to the centre of mass and 
stronger field is experienced by light ray. So we obtain much higher values for the total rotation of the polarization 
vector. We take couple of test cases and let us choose the cases of two millisecond pulsars PSR J 1748-2446ad 
Hessels et al. (2006); Nunez & Nowakowski (2010); Dubey & Sen (2014) and PSR B 1937+21 Ashworth et al. (1983); 
Nunez & Nowakowski (2010); Dubey & Sen (2014) for which the rg and ‘a’ values are (4.050km., 2.42325km.) and 
(4.050km., 2.19438km.) respectively. These two pulsars have been chosen arbitrarily, but we tried our best to choose 
pulsars whose Schwarzschild radii(rg) and rotation parameters (a) values are as close as that of Sun, so that the 
comparison will be easy. We take the impact parameters for them as the physical radii which are 20.1km. Hessels et al. 
(2006) and 20.2km. Ashworth ct al. (1983) respectively. By using Eqn.(35) we find the values of Axeg to be 2°.2471 and 
2°.0172, for the two pulsars PSR J 1748-2446ad and PSR B 1937+21 respectively. So the results are quite appreciable. 
These results have been listed in Table 1. In Fig. 2, we plot Axeq against Xeq for the pulsar PSR J 1748-2446ad for 
various values of impact parameters from the centre of mass. The impact parameter values were taken as 1 time, 1.5 
times and 2 times of the physical radius. We can clearly see the dependence of the rotation of polarization vector Axeq 
n the initial value of polarization vector ( Xeq) ■ We observe this form Fig. 2 rotation value attains a maximum, when 
initial position angle is 47 degrees approximately. From Eqn.(34) also, we get the same value of 47.25 degree for the 
rotation value to be maximum. 

As a further illustration of our calculation, we shall now consider prograde and retrograde orbits of light ray with 
respect to the rotation axis of the gravitational (or Kerr) mass. In case of prograde orbit the equation would be 
same as Eqn.(29). But for the reverse or retrograde orbit, the rotational parameter (a) would be taken as negative in 
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Eqn.(29). Thus the expressions for rotation of polarization vector would be written as: 


and 


Axeqipro) = arctan 



Axeq{retro) = —arctan 




(36) 


(37) 


We note that Eqn.(36)is same as Eqn.(ou) obtained earlier, for prograde. Calculated values for prograde and 
retrograde orbits differ only in sign or direction. 


4.0.3. Case III: Ez = 0 i.e. ^ = oo 

Considering only the Z component is present between two electric field components (xeg=90°). Now going back to 
Eqn.(29) we can show that lim^_>oo Axeq = 0 by applying L’Hospital rule. This is also evident from Eig.(2). So in this 
case there will be no rotation of polarization vector as evident. At this stage, it will be worth to make a compassion 
with the work done by Su & Mallett (1980) earlier. Those authors found that, for an asymptotic flat observer) GIE 
case), the rotation of polarization vector is solely due to the frame dragging effect. This is in perfect agreement with 
our observation in the present work. Eurther, Su Mallett (1980) considered distances beyond 6rg and obtained the 
amount of rotation to be upto few degrees. This is also grossly in agreement with what we have observed. An exact 
comparison between the approaches and results obtained by Su Sz Mallett (1980) with those by us, is not possible due 
to several reasons. The former authors, obtained numerical results by running some specific codes and sometimes a 
very high degree of approximation was applied. Also it appears that, they limited their rotation parameter values (in 
Kerr line element) lower than the Schwarzschild impact parameter. In addition, they obtained solutions under the 
cases where the impact parameter (distance of closest approach of the light ray to the massive body) is at least six 
times the Schwarzschild radius. Eurther, the authors also limited all their calculations, to an initial condition, where 
polarization of light ray is at 45 °with respect to the equatorial plane. 

In our approach, all the above three restrictions are relaxed and the solutions that we obtained are of more general 
in nature. However, the underlying physics in the two approaches appear to be the same. This is obvious from the 
Equation (4.27) of Su & Mallett (1980) and our Eqn.(ll), where the rotation angle (of polarization vector) depends 
on the same set of metric elements and the mathematical nature of the dependence is the same. In spite of the finer 
differences between the two approaches, we can make an attempt to compare the results obtained by the two methods 
as follows. As can be seen from Fig. (3) of Su & Mallett (1980), the rotation angle is about 2.5 degrees ( for GIF) 
under the conditions of rotation parameter= Vg and impact parameter=6*rg . With these same set of parameters, we 
can calculate from our Eqn.(35) the rotation of polarization vector to be about 2.69 degrees. Since our work, has fewer 
assumptions and approximations, we claim the results obtained by us to be more accurate. In addition, the present 
work is more general in nature, direct and has more flexibilities to accommodate different geometries. 


5. POLARIZATION INTRODUCED IN THE LIGHT RAY DUE TO CURVED SPACE-TIME 

In this section, we explore whether the effect of curved space-time introduces any polarization for light ray propa¬ 
gating through it or not? Now from the definition of Stoke’s parameter Stokes (1852); Born & Wolf (1999) for linearly 
polarized electromagnetic wave, one can write the expression for polarization (P) in terms of the Stokes vectors 


Table 1. Schwarzschild radius {rg), Physical radius (R or b), Rotation parameter (a) and Ax^q as the total rotation of 
polarization vector for different millisecond pnlsars and Sun are shown.(original Xeq =45°, C = 1 ) 


Name of Stars 

Tg (km.) 

R(km.) 

a (km.) 

AXeq 

Sun 

3.000 

0.7x10® 

1.69749 

4".2873 X 10"® 

PSR J 1748-2446ad Hessels et al. (2006) 

4.050 

20.1 

2.42325 

2°.2471 

PSR B 1937-1-21 Ashworth et al. (1983) 

4.050 

20.2 

2.19438 

2°.0172 
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(StokesStokes (1852)) as: 

I=<Dl> + 
Q=< Dl> - 
U = 2 < DyD, 


<Dl> 

(38a) 

<Dl> 

(38b) 

cos <5 > 

(38c) 


where S is the phase difference between Dy and and < .. > indicates the time average. In the present work, it is 
clear that that both the components Dy and are in same phase. From Eqn.(10a)and Eqn.(lOb) we can calculate 
Dy and Dz, where it has been assumed that the light is linearly polarized. So from the above definition of polarization, 
we can write(omitting < .. > sign everywhere for our convenience). 


^{Dl-Dir + {2DyDzr 
I DI + Dl 


(39) 


Thus we find that, the light remains same way 100% plane polarized, as it was before entering the gravitational field. 
As a result, one can conclude here that the curved space time geometry (of Kerr field), does not introduce any new 
polarization in the electromagnetic wave (or light). However, the gravitational field just rotates the direction of the 
polarization vector (or electric vector). 


6. RESULT AND DISCUSSIONS 

The above work clearly demonstrated and systematically calculated, the effect of curved space-time on the polariza¬ 
tion states of a light ray. From this work, one can easily come to a conclusion that there is no effect on polarization, 
due to the gravitational field of static non-rotating body. For example the Schwarzschild field will have no effect on 
the polarization of the electromagnetic field so far, which is very obvious and supports the previous work by several 
authors. 

Present work clearly indicates that, however the gravitational field generated by a rotating body would affect the 
direction of polarization vector of electromagnetic wave. The work presented here provides an analytical expression to 
find out the amount of rotation of polarization vector. This expression is much simpler and straightforward compared 
to the ones obtained by previous authors ; viz the work by Su & Mallott (1980) requires many approximations and 
numerical techniques to find the exact value of rotation. Calculations carried out here, for couple of pulsars as test 
cases, demonstrated that the result is appreciable and it is in the order of few arc degrees, which can be measured by 
present day telescopes. However, for Sun the calculated rotation is in the order of few micro-arcseconds, and this can 
be measured by new generation telescopes and polarimeters. Also based on the calculations presented here, conclusions 
can be drawn that the gravitational field would not produce or introduce any polarization in the light and it will only 
rotates the direction of existing polarization vector. 

Our present study may also open new avenues for the indirect detection of gravitational waves (Abbott et al. (2016)). 
This may be possible in two different ways (i) the electromagnetic counterpart of gravitational wave may be emitted in 
a way involving axial symmetry [viz. binary black hole coalescence Abbott et al. (2016)), which will cause the rotation 
of polarization vector of electromagnetic wave passing close by; and/or (ii) the electromagnetic wave passing through a 
region of gravitational wave will be lensed Sbytov (1973), resulting in a good possibility of rotation of its polarization 
vector, especially if the gravitational wave is plane polarized. However, these are only possibilities, which we would 
like to explore in our future work. 
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